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ABSTRACT 

Some massive binaries should contain energetic pulsars which inject relativistic lcptons 
from their inner magnetospheres and/or pulsar wind regions. If the binary system is 
compact enough, then these leptons can initiate inverse Compton (IC) e ± pair cascades 
in the anisotropic radiation field of a massive star. 7-rays can be produced in the IC 
cascade during its development in a pulsar wind region and above a shock in a massive 
star wind region where the propagation of leptons is determined by the structure of a 
magnetic field around the massive star. For a binary system with specific parametres, 
we calculate phase dependent spectra and fluxes of 7-rays escaping as a function of the 
inclination angle of the system and for different assumptions on injection conditions of 
the primary leptons (their initial spectra and location of the shock inside the binary) . 
We conclude that the features of 7-ray emission from such massive binaries containing 
energetic pulsars should allow to obtain important information on the acceleration of 
particles by the pulsars, and on interactions of a compact object with the massive star 
wind. Predicted 7-ray light curves and spectra in the GeV and TeV energy ranges from 
such binary systems within our Galaxy and Magellanic Clouds should be observed by 
future AGILE and GLAST satellites and low threshold Cherenkov telescopes such as 
MAGIC, HESS, VERITAS or CANGAROO III. 

Key words: gamma-rays: massive binaries - gamma-rays: theory - radiation mech- 
anisms: nonthermal 



1 INTRODUCTION 



It seems obvious that some close massive binary systems 
should contain young neutron stars able to accelerate lep- 
tons and possibly ions to relativistic energies. In fact, nu- 
merical simulations of the evolution of neutron stars show 
that the fraction of systems with not accreting pulsars (dur- 
ing so called ejector phase) may become several percent of 
the total number of massive binaries (Lipunov 1990). How- 
ever, due to their proximity to the high density winds of 
the massive stars most of such systems do not show clear 
modulation of their radio signal with pulsar rotational pe- 
riod. Therefore, only in the case of broad or highly eccentric 
binaries the pulsars with periods significantly shorter than 
1 second have been discovered (e.g. PSR B 1259-63 with the 
period of 47.8 ms or A0538-66 with the period of 69.2 ms). 
TeV 7-ray emission, modulated with the period of ~ 12.59 
ms, has been also claimed from the compact binary Cyg X-3 
(Brazier et al. 1990), but this result has not been confirmed 
by any other experiment. Other binaries are suspected to 
contain short period pulsars (e.g. LSI +61 303, Maraschi & 
Treves 1981), or stellar mass black holes ejecting relativistic 



particles highly anisotropically (e.g. Cyg X-l, Bednarek et 
al. 1990). 

In fact, observations of high energy X-ray and 7-ray 
sources to be coincident with the locations of some binary 
systems support the hypothesis that the high energy pro- 
cesses play an important role in these sources. For example, 
EGRET sources have been found in directions towards e.g. 
LSI +61 303 (2EG J0241+6119, Thompson et al. 1995), Cyg 
X-3 (2EG J2033-4112, Mori et al. 1997), and LS 5039 (3EG 
J1824-1514, Paredes et al. 2000). However, no TeV 7-rays 
from the above mentioned sources have been observed up to 
now (except early claims during the 80's which were not 
confirmed by subsequent more sensitive observations; see 
e.g. Weekes 1992). The limits on TeV emission from PSR 
B1259-63, recently observed by the CANGAROO group at 
~47 days and at ~157 days after the periastron (Kawachi 
et al. 2004), are above the theoretical predictions (Kirk et 
al. 1999, Kawachi et al. 2004, Murata et al. 2003). However, 
very recently the HESS group reported positive detection of 
this binary system on the level of 5% of the Crab emission 
at energies above ~ 400 GeV with the power law spectrum 
and spectral index -2.8 (Schlenker et al. 2004). The TeV 
emission has been observed a few days before and after the 
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periastron passage and declines towards the periastron mo- 
ment. Observations by the Whipple telescope of other four 
binary systems containing young pulsars (Hall et al. 2003) 
have not shown any positive detection of a steady or modu- 
lated TeV signal. Another binary system, Cen X-3, reported 
more recently in the TeV 7-rays (Chadwick et al. 1998, 1999, 
Atoyan et al. 2002), and also in the GeV 7-rays (Vestrand 
et al. 1997), contains slowly rotating and accreting neutron 
star. 

In spite of not completely convincing high energy obser- 
vational results, massive binaries are still often considered 
as likely sources of 7-rays produced mainly in two general 
scenarios: (1) the anisotropic injection of particles from neu- 
tron stars or black holes, which interact with the radiation 
or matter inside the binary system (a massive companion, 
an accretion disk and its corona or a blob, see e.g. Cheng 
& Ruderman 1991, Aharonian & Atoyan 1991, Levinson 
& Blandford 1996, Aharonian & Atoyan 1999, Romero et 
al. 2001, Atoyan et al. 2002, Georganopoulos et al. 2002, 
Romero et al. 2002, Romero et al. 2003, Bosch-Ramon & 
Paredes 2004a,b, Orellana & Romero 2004), or (2) the inter- 
action of particles accelerated by a compact object or a shock 
wave created in collisions of the pulsar and stellar winds 
(e.g. Vestrand & Eichler 1982, Harding & Gaisser 1990, Ta- 
vani et al. 1994, Tavani & Arons 1997, Kirk et al. 1999, 
Ball & Kirk 2000, Murata et al. 2003) or two stellar winds 
(e.g. Eichler & Usov 1993, Benaglia & Romero 2003). In 
fact, the injection of relativistic particles and high energy 
7-rays in dense radiation field expected close to the ac- 
cretion disks around compact objects (Carraminana 1992, 
Bednarek 1993) and massive companion stars (Protheroe & 
Stanev 1987, Moskalenko et al. 1993) should result in co- 
pious 7-ray production. In very compact binaries and com- 
pact objects surrounded by luminous accretion disks, the 
primary particles can initiate IC e ± pair cascades triggered 
due to the large optical depths. The 7-ray spectra are then 
produced with the characteristic cut-offs at a few tens of 
GeV (Bednarek 1997, Bednarek 2000). It was assumed in 
the latter works that the secondary e pairs produced in 
the cascade process are isotropised inside the binary system 
by the random component of the magnetic field. 

In the present paper we consider in detail a more spe- 
cific scenario in which the volume of the binary system is 
separated by the shock wave into two regions with differ- 
ent properties. The shock appears as a result of collisions of 
the pulsar and massive star winds. We analyze very com- 
pact binary systems which contain a young pulsar and a 
massive companion star of the OB or WR type able to cre- 
ate soft radiation field in which the optical depths for rel- 
ativistic leptons are much larger than unity. In fact, young 
pulsars in binary systems can be responsible for quite dif- 
ferent phenomena such as the appearance of jets called mi- 
croquasars (e.g Cyg X-3), binary radio pulsars (e.g. PSR 
1259-63, SAX J0635+0533) or a long period accreting ob- 
jects (e.g. A0535+26). It is very difficult to observe some of 
such systems due to dense stellar winds produced by early 
type massive stars. However, binaries hidden dense winds 
should manifest themselves by the presence of large non- 
thermal luminosities which can, as we argue below, peak in 
the 7-ray energies. In Sect. 2 we define general conditions 
inside the massive binary in which the termination shock is 
created in collisions of the pulsar and stellar winds. The in- 



jection of relativistic leptons (electrons or e pairs) by the 
pulsar, their subsequent propagation and interaction with 
the radiation of the massive star is considered. In Sect. 3 we 
calculate production of 7-rays in IC e ± cascades initiated by 
these leptons in the pulsar wind zone (PWZ) by applying 
the Monte Carlo method. Because of partial absorption of 
these 7-rays in the massive star wind region (MSWR), the 
secondary e ± pairs are created and initiate there the second 
part of the cascade. This cascade differ significantly from the 
cascade in the PWZ due to the presence of ordered magnetic 
fields. The spectral and angular features of 7-rays escaping 
from the PWZ and MSWR are discussed in Sect. 5. As an 
example we consider the binary system with the parameters 
derived for the WR star in the Cyg X-3 binary in which 
there are some observational evidence of the existence of a 
very young and energetic pulsar with the period of several 
milliseconds. 



2 A PULSAR CLOSE TO A MASSIVE STAR 

We consider the binary system containing a fast pulsar and 
a massive companion star (M « 10M Q ). It is assumed that 
the energy loss rate of the pulsar, L rot , is high enough that 
the matter from the massive companion can not accrete 
onto the pulsar surface neither from the outflow through 
the Lagrangian point or from the dense stellar wind. This 
means that the binary system has to be in the ejector phase, 
in which the structure of the inner pulsar magnetosphere, 
where efficient acceleration and e bfpair creation occur, is 
not destroyed. The binary systems having the pulsar with 
the period fulfilling the condition (Harding & Gaisser 1990), 

P ms < 3lBti 7 Mr 8 2/7 , (1) 

should blow away accreting material, where P ms is the pulsar 
period in ms, B = 10 12 £>i2 G is its surface magnetic field, 
and M = 10 18 Mis g s _1 is the accretion rate. 

Let us assume that the pulsar is on a circular orbit 
around the massive star with the radius R s , effective surface 
temperature T B , and surface magnetic field B s . The star cre- 
ates the wind which has the termination velocity Vao and 
is characterized by the mass loss rate M. The separation 
of the stars is D. As a result of the interaction of the pul- 
sar and stellar winds a double shock structure is formed, 
separated by the contact discontinuity, at the distance de- 
termined by the above mentioned parameters of the stars. 
We apply the simplified model for the structure of the col- 
liding winds based on momentum conservation (Girard & 
Wilson 1987). In this model the shock structure reaches a 
steady state configuration under the assumption of negligi- 
ble pulsar orbital motion. In this case, the shock geometry 
of specific binary system is determined by the quantity 77 
(Ball & Dodd, 2000), defined as follows: 

n = L ro t/(MV w c), (2) 

where the wind velocity, V w , depends on the radial distance, 
r, from the massive star according to (Hamann 1985) 

Vw(r) =« 00 (l-P A /r) 1 ' 5 , (3) 

where Ra is the Alfven radius which can be derived by solv- 
ing the equation (1 — R b /Ra) = £,(R b /Ra) 4 - It has simple 
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Figure 1. The schematic scenario of the interacting neutron star 
and massive companion inside the compact binary system. Pri- 
mary leptons, with the Lorentz factors f e , are injected by the 
pulsar and propagate inside the pulsar wind region (PWZ) along 
direction defined by the angle 0, comptonizing soft radiation from 
the massive star. Secondary 7-rays can be absorbed in the same 
radiation field either inside the PWZ or, after passing the pulsar 
wind termination shock, in the massive star wind region (MSWR), 
triggering an electromagnetic cascade. Leptons propagate radially 
inside the PWZ but follow the magnetic field structure inside the 
MSWR. 7-rays from the cascade either escape from the binary 
system or fall onto the surface of the massive star, depending 
on the injection parameters of the primary leptons (their ener- 
gies and directions) and the parameters of the binary system (the 
massive star parameters and system separation). 



approximate solution 



Ra = R B X 



1 e /4 , 



e» 1, 



(4) 



where £ = B a 2 : R S / '(Mvoo)- 

The closest radial distance from the pulsar to the termina- 
tion shock (measured in the plane of the binary system) is 
then equal to, 



Po = D y/Tj/{l + y/Tj). 



(5) 



In order to estimate po we have to solve the set of Eqs. 2-5. 
Note, that for rj < 1 the star wind dominates over the pulsar 
wind and the termination shock wraps around the pulsar. 
For r) > 1 the shock wraps around the companion star. For 
large distances from the massive star the shock tends asymp- 
totically to a cone characterized by a half-opening angle tp, 

v> = 2.1(1 -7r 2/5 /4)»r 1/3 . (6) 

We approximate the structure of the termination shock by 
the sphere of radius po for 8 > 80 (0 is defined in Fig. 1), and 
by a cone for < 8 < 9o, where the angle do is defined as 
8 and is given by the perpendicular line from pulsar to the 
generator of the cone. The shock is unterminated for 8 < tp. 

It is assumed in these simple estimations that the pulsar 
and stellar winds are symmetric. In a more realistic case the 
dependence of the wind pressure on the distance from the 
equatorial plane of both stars should be taken into account. 
This might result in more complicated structure of the shock 
surface and involve more parameters which would prevent 
clear analysis of the angular features of high energy 7-ray 



emission from such compact binaries. In a future work (Sier- 
powska & Bednarek 2004a), in which we apply such a model 
to the 7-ray production in the binary pulsar PSR 1259-63 
on the orbit around the Be star, the effects associated with 
the non-spherical stellar wind will be taken into account. 

In order to perform detailed calculations in the case 
of specific binary we apply as an example the parameters 
expected for the massive WR star in the Cyg X-3 system. 
This is a short period compact binary, r = 4.8 hr, with 
the massive star radius < 3 — 6Rq, surface temperature 
T s > (7 — 9) x 10 4 K, separation of the components 3.2 < 
D/Rq < 5.6, and the mass loss rate M ~ l.lxlO _5 M Q yr' 1 
(Cherepashchuk & Moffat 1994). The orbital, phased re- 
solved infrared spectra in the outburst and quiescent stages 
are consistent with the high orbital inclination of the sys- 
tem with respect to the observer at the Earth, i > 60°, 
provided that the mass of the massive star is in the range, 

5 < Mwr/M q < 11 (Hanson et al. 2000). The infrared ob- 
servations by van Kerkwijk et al.(2002) are consistent with 
the high orbital inclination, i = 74°, and the mass loss rate 
in the range M ~ 1.2 x 10~ 4 M© yr _1 (based on the infrared 
data) and M ~ 0.6 x 1O~ 5 M0 yr _1 (based on the increase 
of orbital period). The stellar wind has terminal velocity 
Voc = 1-45 x 10 3 km s _1 . On the other hand, the inclina- 
tion angle of system derived from the X-ray Chandra data 
is much lower, i « 24°. The mass of the compact object has 
been constrained by < 3.6M Q , and the mass and radius of 
the stellar companion by < 7.3M© and i? s < 1.6i?o (Stark 

6 Saia 2003). 

Based on the above observational constraints, we adopt 
the following parameters for our specific binary system. A 
Wolf-Rayet type star with a radius R B = 1.6 x Rq, an effec- 
tive temperature, T e g = 1.36 x 10 J K, and a typical surface 
magnetic field, B s ~ 10 2 -10 3 G (Usov & Melrose). The mass 
loss rate is M ~ 0.8 - 8.0 x 10~ 5 M yr -1 , the velocity of 
the stellar wind at infinity v x ~ (1 — 5) x 10 8 cm s~\ and 
the star rotational velocity v ro t ~ (0.1 — 0.2)voo. The pul- 
sar has a period P ms = 12.59 and a surface magnetic field 
B s = 4.95 x 10 11 G (Brazier et al. 1990). It is on a circular 
orbit, with a separation of D = (3.6 ± 1.2)7?© = 2.25 x R B 
(Stark & Saia 2003). For these parameters, the pulsar energy 
loss rate is, 



6 x 10 43 Bl 2 



erg 



6 x 10 38 erg s" 



(7) 



For the above mentioned parameters of the massive compan- 
ion in the Cyg X-3, the value of the parameter r\ (see Eq. 2) 
is in the range 0.067 < jy < 0.67 (for fixed V x = 1.45 x 10 3 
km s _1 ). For further calculations we apply the average value 
rj = 0.3 and for comparison the lowest value r\ — 0.06. 

The conditions in the two regions of the binary system, 
i.e the PWZ and the MSWR, separated by the shock struc- 
ture differ significantly. Below the shock, inside the PWZ, 
relativistic leptons are frozen in the magnetized pulsar wind 
which propagate radially from the pulsar. Therefore, their 
synchrotron losses are neglected in our cascade calculations. 
We consider the IC cascade in the PWZ assuming that the 
thermal radiation from the massive star dominates in the 
main volume of this region, neglecting the other sources of 
soft photons e.g. such as the thermal emission from the neu- 
tron star surface or the nonthermal emission from the in- 
ner pulsar magnetoshere. This is true if the temperature of 
the whole neutron star surface is below Tns ~ 6 x 10 6 K. 
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But such surface temperature is characteristic for the cool- 
ing neutron star in the age of only a few months after for- 
mation (Nomoto & Tsuruta 1987). This limit temperature 
has been obtained by comparing the energy density of ther- 
mal radiation from the massive star, U B w assT* (R B / D) 2 , 
where <tsb is the Stefan-Boltzman constant, with the en- 
ergy density of thermal radiation from the neutron star, 
(7ns ~ o"sb2n S (_Rns/-Rlc) 2 (1 + cosfl), at the distance of 
the light cylinder radius, _Rlc- For the period of the pul- 
sar as suggested in the Cyg X-3 binary, P — 12.59 ms, we 
obtain Rns/Rlc ~ 60. The approximate angle of interac- 
tion between leptons and photons coming from the neutron 
star surface is then 9 = n — Rns/Rlc radians, assuming the 
case of a spherical wind. The estimate of the possible con- 
tribution of the nonthermal radiation from the pulsar inner 
magnetosphere is very difficult since the geometry of this 
nonthermal emission is not well known. We comment that 
this nonthermal photons has to be strongly collimated along 
the direction of motion of relativistic leptons as postulated 
by the observation of very narrow peaks in the light curves of 
young pulsars. Therefore, its importance should be strongly 
suppressed. However, Bogovalov & Aharonian (2000) calcu- 
lated possible 7-ray emission from the vicinity of the light 
cylinder in the case of the Crab pulsar applying a specific 
model for the radiation close to the light cylinder radius. Ac- 
cording to these calculations the optical depth for leptons in 
the thermal radiation field of the Crab pulsar is very low 
at the light cylinder. However, in the nonthermal radiation 
from the pulsar inner magnetosphere, the optical depth for 
leptons can be significant at the distance up to a few light 
cylinder radii. We conclude that the soft radiation in the 
main part of the PWZ is dominated by the thermal radiation 
from the massive star. Although, the important contribution 
of the nonthermal radiation from the inner pulsar magneto- 
sphere at distances just above the light cylinder radius can 
not be excluded. This problem should be studied in a more 
detail when reliable models of the high-energy processes in 
the pulsar magnetospheres are better constrained. 



Leptons which move through the PWZ interact effi- 
ciently with the soft radiation of the massive companion 
initiating IC e pair cascade. The charged products of this 
cascade arrive finally to the shock region in the pulsar wind 
and follow the flow along the shock surface. The power in 
these secondary leptons is relatively low with respect to the 
power in secondary cascade 7-rays in the case of very close 
binary systems considered in this paper. The secondary 7- 
rays move into the massive star wind region. Some of them 
escape the binary system but a significant part can be con- 
verted into the next generation of e ± pairs which have to 
follow the complex structure of the magnetic field present in 
the stellar wind. These pairs can trigger further cascading 
processes producing next generation of 7-rays at directions 
which depend not only on the injection geometry of primary 
leptons but also on the geometry of the magnetic field. All 
these processes are discussed in detail in this paper and the 
spectra of 7-rays escaping at different angles with respect to 
the plane of the system and as a function of the phase of the 
pulsar on its orbit around the massive star are calculated. 



2.1 Injection of leptons by the pulsar 

Young pulsars are efficient sources of energetic leptons, i.e. 
electrons and positrons. These leptons originate in the cas- 
cade processes occurring in the inner pulsar magnetosphere, 
as expected in terms of the polar cap model (e.g. Ruderman 
& Sutherland 1975, Arons & Sharlemann 1979, Daugherty 
& Harding 1982) and the outer gap model (e.g. Cheng et 
al. 1986). Leptons escape through the pulsar light cylinder 
to the pulsar wind zone (PWZ) where they can be addi- 
tionally accelerated either just above _Rlc (e.g. Beskin & 
Rafikov 2000) or linearly through the PWZ (e.g. Contopou- 
los & Kazanas 2002), or at the pulsar wind shock (e.g. Ken- 
nel & Coronity 1984). 

In this paper we consider two models for the primary 
spectra of leptons injected into the radiation field of the 
massive star: 

(i) The power law spectrum of injected leptons between 
100 MeV and 500 GeV and the spectral index -1.2, as en- 
visaged in the recent calculations of the spectra of leptons 
escaping from the inner magnetopshere performed by Hib- 
schman & Arons (2001). 

(ii) The monoenergetic injection of leptons with energies 
10 6 MeV corresponding to the Lorentz factors of the pulsar 
wind with the parameters typical for the Crab pulsar. These 
leptons have very similar energies to those ones expected 
for the supposed pulsar in Cyg X-3 binary system due to 
the similar value of B/P 2 , which determine the maximum 
potential drops through the pulsar magnetopshere. In fact, 
adopting the formulas of the slot-gap pair-cascade model 
(Arons 1983) for a e pair generation rate, N e ± , and small 
ratio <7pw of Poynting flux to kinetic energy flux in the wind 
not far from the light cylinder radius, the Lorentz factor 
7 pw of the pulsar wind can be estimated from (Melatos et 
al. 1995), 

7 P w = Lrot/(l +cr pw )iV e ±m e ±c 2 , (8) 

where m e ± is the electron mass, and c is the velocity of light. 
For the considered here parameters of the pulsar N e ± ~ 
2.8 x 10 37 s _1 and the pulsar wind Lorentz factor can be as 
high as 7 pw ~ 10 7 (Melatos et al. 1995). The small ratio a pw 
can be obtained in the narrow boundary layer near the light 
cylinder as suggested by Beskin & Rafikov (2000). These 
authors argue that the energy of the magnetized wind may 
be transfered to leptons with the Lorentz factors ~ 10 6 . 

Most of the calculations presented below are shown for 
the power in injected spectrum of leptons normalized to 10 6 
MeV sr _1 in order to allow direct comparison between the 
cases of power law and monoenergetic injection of leptons. 
The absolute fluxes of 7-rays produced in the cascade can 
be then obtained by multiplying by the factor ~ 3 x 10 37 s _1 
which is calculated from normalization of the power equal 
to 10 6 MeV to the total rotational energy lost by the pulsar 
(see Eq. 8). 

2.2 Conditions for propagation of leptons and 
gamma-rays 

We are interested in binary systems in which relativistic 
leptons, injected by the pulsar, can develop the IC e^ 1 pair 
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Figure 3. The 7-ray spectra produced in the cascade inside the PWZ in the case of injection of primary leptons with: the monoenergctic 
spectrum and energy 10 6 MeV (on the left), and the power law spectrum (on the right). The spectra are shown for different injection 
angles of the primary leptons: = 0° (thick solid curve), 30° (dashed), 60° (dot-dashed), 90° (dotted), 120° (thin solid), and 150° 
(dot-dot-dot-dashed). The 7-ray spectra which arrive to the shock region (defined by r\ = 0.3) are shown in (a) and (c), and the 7-ray 
spectra escaping to the observer at the infinity (after their absorption inside the MSWR) are shown in (b) and (d). The spectra are 
obtained by sorting photons within the intervals with the width A (log E) = 0.2 with the lower boundary of the interval marked in the 
figures. 



cascades in the anisotropic radiation of the massive compan- 
ion. Considered by us binary system meets this requirement 
if the optical depths for IC scattering of massive star thermal 
radiation and for 7-rays absorption process in the same ra- 
diation are much larger than unity. For the injection place of 
leptons and 7-rays equal to the separation distance between 
the pulsar and the WR star, D — 2.25R S , and other param- 
eters as mentioned above, we calculate the optical depths up 
to the infinity for different angles of injection measured with 
respect to the direction defined by the centres of these stars 
and as a function of lepton energies (see Fig. 2). The optical 
depths for leptons are very large in considered energy range. 
Moreover, 7-rays with energies above a few GeV have high 
probability of interaction for all directions (optical depth 
larger than unity). Therefore, leptons injected with energies 
large enough should develop IC e ± pair cascade. 

In fact, the cascade can start to develop already inside 
the PWZ (see Fig. 2b). The optical depths for all injection 
angles, calculated up to the location of the shock defined by 
rj = 0.3, are much larger than unity at specific energy ranges. 
In Fig. 2c the dependence of the optical depths for leptons 
and 7-rays on the location of the shock inside the binary 
system is shown. It is clear that even for the extreme values 
of 77 = 0.06 and 0.6 leptons should develop efficient cascades 
if injected with energies above ~10 GeV. ^,From comparison 
of Figs. 2a and 2bc we find that the optical depths for 7-ray 
photons which are produced below the shock (in the PWZ) 
but propagate through the MSWR are still much larger than 
unity. Therefore, these 7-rays partially convert into the next 
generation of e pairs in the MSWR. These e^ 1 pairs pro- 
duce next generation of 7-rays which angular distribution on 
the sky may differ from initial injection directions of primary 



leptons due to the presence of a complex magnetic field in 
the MSWR. For that reason we consider the cascades devel- 
oping inside the PWZ and MSWR separately. 



3 PULSAR WIND REGION 

It is assumed that energetic leptons are injected from the 
vicinity of the pulsar light cylinder and propagate radially 
from the pulsar almost at rest with respect to the pulsar 
wind. In this case we neglect the synchrotron losses of lep- 
tons during their propagation up to the pulsar wind shock. 
If the optical depths are such as shown in Figs. 1, leptons 
develop IC e ± pair cascade whose efficiency depends on the 
injection parameters of these primary particles (initial ener- 
gies of leptons, their injection directions and parameters of 
the binary system and massive star). We assume that the 
cascade initiated by specific primary lepton in the PWZ de- 
velops in one dimension, i.e. in the direction of propagation 
of the primary particle. Such one-dimensional cascade de- 
velops up to the pulsar wind shock. The secondary 7-rays 
pass through the shock into the massive star wind region 
where they can be additionally absorbed producing next 
generation of energetic e ± pairs. For considered large opti- 
cal depths, only the 7-rays with energies below the threshold 
for e pair production escape freely from the binary system. 
The secondary e ± pairs which are produced in the PWZ are 
captured by the magnetic field of the shock region and move 
along its surface with the pulsar wind plasma. However, in 
the case of compact binary systems such as Cyg X-3, we can 
neglect the contribution of secondary e pairs to the total 
escaping 7-ray spectrum, because the part of energy of the 
primary leptons transfered to the secondary e^ 1 pairs is rela- 
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Figure 2. The optical depths for leptons on ICS process, and 
for 7-rays on eP 1 pair production in the anisotropic radiation of 
the massive star during their rectilinear propagation from the 
injection place at the distance D = 2.25R S from the massive star 
up to the infinity (figure (a)). The optical depths are shown as a 
function of particle energies for the case of injection at the angle 9, 
measured from direction defined by the centers of these stars (see 
Fig. 1): = 0° (thick full curves), 30° (dashed), 60° (dot-dashed), 
90° (dotted), 120° (thin full), and 150° ( dot-dot-dot-dashed), (b) 
As in (a) but for particles propagating only to the termination 
shock defined by r\ = 0.3 and for the injection angles 90° (solid 
lines), 120° (dotted lines) and 150° (dashed lines). Note that for 
the angles 0° — 60°, the optical depths are the same as in (a) due 
to the lack of boundary on the rectilinear propagation caused 
by the presence of the shock, (c) As in (b) but for two specific 
locations of the shock defined by r) = 0.06 and 0.6 and for the 
injection angles of primary leptons equal to 90° (solid lines) and 
150° (dashed lines). 

tively low compared to that one transfered to the secondary 
7-rays, as we show later. The influence of this part of lep- 
tonic cascade on the total 7-ray spectrum escaping from the 
system will be considered in the future work dedicated to 
wider binary systems. 

3.1 Gamma-ray spectra from the PWZ 

We follow the development of such IC e^ 1 pair cascade in 
the anisotropic radiation of the massive star assuming that 



the primary particles are injected in the place corresponding 
to the location of the pulsar in the binary system (the light 
cylinder radius can be safely neglected with respect to the 
characteristic dimension of the PWZ). The procedure for the 
cascade Monte Carlo simulations in the PWZ is generally 
the same as used in Bednarek (1997, 2000) except for the 
assumption on the local isotropisation of secondary e^ 1 pairs 
applied in that paper. In this paper it is assumed that the 
secondary cascade e^ 1 pairs follow the direction of the parent 
7-rays. The distance of the first interaction of a lepton in 
the IC scattering process, Lie, is determined by sampling 
from the optical depth calculated for specific parameters of 
lepton injection place. Lie is found for the random number 
Pi applying the formula, 

Pi = exp (- J IC Aic 1 (E, L, 9)dl)j , (9) 

where Aic(-E, L, 9) is the mean free path for IC scattering of 
lepton with energy E, at the place defined by the propaga- 
tion distance L and the angle 6 (see Fig. 1). Aic is calculated 
from the formula 

\Tc(E,L,9) = J dn{l + n) J d<f> j N(e,Q)a K Nde, (10) 

where fi — cos 8 is the cosine of the angle between lepton and 
the incident soft photon, N(e, fi) is the differential density of 
soft photons, and <tkn is the Klein-Nishina cross section. The 
limits of integrations in this formula are constrained by the 
geometry of the scattering process (propagation directions 
of leptons and thermal photons from the massive star) and 
the kinematics of the IC process. Since these limits are very 
complicated (several cases have to be considered), they will 
be published separately in Sierpowska (2004). The distance 
at which secondary e ± pair is produced in collisions of the 
7-ray photon with the soft photon is obtained in a similar 
way (see Eqs. 9 and 10) by replacing the Klein-Nishina cross 
section by the cross section for pair production in the two 
photon collision (see Appendix A in Bednarek 1997). The 
subsequent places of the interaction of secondary leptons 
and 7-rays are determined in the same way by replacing the 
lower limit in the integral in Eq. (9) from the n-th cascade 
step by the upper limit in the integral from the previous, (n- 
l)-th cascade step. The neutron star is relatively close to the 
massive star with respect to its dimensions as the separation 
of the components is only 2.25 larger than the radius of the 
massive star. Therefore, the geometrical effects connected 
with limited solid angle of impinging soft photons have to 
be taken into account when considering the IC scattering 
process. 

The energies of secondary particles produced by par- 
ent leptons or 7-rays as a result of the interaction with the 
soft photons from the massive star are obtained for the ran- 
dom number P2 by sampling from the differential energy 
distributions of these secondary particles according to the 
relation 

P 2 = ^ NdE^j ^ NdE^j , (11) 
where N is the differential spectrum of secondary particles 
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30 
60 
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(b) 
30 
60 
90 

120 90 74 84 52 

150 87 64 85 56 

180 87 - 85 - 



Table 1. The percentage of primary lcptons initial energy with 
monoenergetic (a) and power law spectra (b) (see text for details) 
which is transfcred to the 7-rays in the cascade process for two 
locations of the shock defined by r\ = 0.3 and 0.06. The part of 
energy in 7-rays which reached the shock region is marked by 
shock and which menaged to escape to the infinity is marked by 
esc. For small angles 9 the shock is not present, and for 8 = 180° 
7-rays do not escape but collide with the surface of the massive 
star. 

{e^ 1 pairs or 7-rays), produced by the parent particle with 
energy E(7~ray or lepton, respectively), calculated at a spe- 
cific location inside the PWZ which is defined by Lie for the 
IC process and by L 7 _ 7 for 7-ray absorption process (see 
for details the Appendix B and C in Bednarek 1997); -Esec.m 
is the maximum energy of produced secondary particle al- 
lowed by the kinematics of specific process, i.e. IC scattering 
or 7 — 7 absorption; E BCC is the simulated energy of produced 
secondary particle. The energy of parent lepton after n-th 
interaction and the energy of secondary leptons, produced in 
7 — 7 absorption, is found by applying energy conservation, 
i.e E e ± t7l — E e ±^ n _ 1 — E~. hn and -E 7 ,„_i = E e + t „ + E e - n . 
Since we are interested in high energy 7-rays, i.e > lQQMeV , 
the cascading procedure is switched off when the secondary 
leptons are cooled to energies below 500 MeV. Such leptons 
are not able to produce 7-rays above 100 MeV in the radia- 
tion described by the black body spectrum with temperature 
~ 10 s K, typical for the massive stars considered by us. The 
simulation procedure of the 7-ray photon interaction place 
and the energies of secondary leptons from photon absorp- 
tion is analogous to what was described above. 



The cascade processes occurring inside the PWZ are 
axially symmetric with respect to the direction defined by 
the centers of stars. Therefore, the spectra of 7-rays, these 
at the location of the shock front and those ones escaping 
to the infinity through the radiation of the MSWR, are cal- 
culated only as a function of the injection angle 9. The 7- 
ray spectra have been obtained for two, mentioned in Sect. 
2.2, distributions of primary leptons and the location of the 
shock defined by r) — 0.3 (see Figs. 3 for the monoener- 
getic and the power law spectra of primary leptons). Note 
that photons are sorted within the intervals with the width 



<t> [°] <t> [°] 

Figure 4. The 7-ray light curves in two energy ranges, 0.1 — 10 
GeV ((a) and (c)) and 10 - 10 3 GcV ((b) and (d)), produced 
in the cascade process occurring inside the PWZ. The cascade 
is initiated by leptons with the monoenergetic (figures on the 
left) and the power law spectra (on the right) for selected values 
of the inclination angles i = 30° (solid curve), 60° (dashed), and 
90° (dotted) measured from the normal to the plane of the binary 
system. 



A(logE) — 0.2. The basic features of the 7-ray spectra ob- 
served at infinity, i.e. after their partial absorption in the ra- 
diation field of the MSWR, can be easily understood if ana- 
lyzed in the context of the optical depths shown in Figs. lab. 
Only primary leptons propagating in the outward direction 
(with respect to the massive star) can produce 7-ray fluxes 
above ~ 100 GeV, being potentially detectable by telescopes 
operating at very high energies (VHE) (Figs. 3bd). On the 
other hand, the 7-ray spectra escaping in the inward direc- 
tions, i.e. close to the limb of the massive star, have compa- 
rable intensities at energies below 10 GeV (HE range) within 
a factor of 2-3. Moreover, these spectra show cut-offs close to 
~ 10 GeV which are determined by the surface temperature 
of the massive star. Therefore, the anticorrelation between 
the 7-ray emission in the HE and VHE ranges is expected 
from such compact luminous binaries. The spectra of 7-rays 
produced in the PWZ and arriving to the shock region are 
shown in Figs. 3ac. The difference between 7-ray spectra 
shown in Figs. 3ac and 3bd gives us the information about 
the importance of absorption processes in the MSWR. It is 
evident that the high energy part of the 7-ray spectrum pro- 
duced in the PWZ is strongly absorbed during propagation 
through the MSWR if the primary leptons are injected in the 
hemisphere containing the massive star, i.e. for the angles 
6 > 90°. Important information on the efficiency of 7-ray 
production can be derived from analysis of the percentage 
of the energy of primary leptons which is converted in the 
PWZ cascade process into the 7-ray photons. Table 1 reports 
these efficiencies at the infinity (marked by esc) and at the 
location of the shock inside the binary system (marked by 
shock) defined by r\ — 0.3 and 0.06. Most of the initial lepton 
energy can be converted into the 7-rays in such a cascade 
already at the distance of the shock in the case of dense soft 
radiation field of such compact massive companion. 

A significant part of the power of produced 7-ray for 
the directions of primary leptons towards the massive star 
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log(E / MeV) log(E / MeV) 

Figure 5. The 7-ray spectra escaping from the binary system for different phases of the pulsar on its orbit around the massive star 
(measured with respect to the location of the observer), ip = 30° (solid curve) 60° (dashed), 90° (dot-dashed), 120° (dotted), 180° 
(dot-dot-dot-dashed), and different inclination angles of the binary system i = 30° ((a) and (d)), 60° ((b) and (e)), and 90° ((c) and (f)). 
7-rays are produced in the PWZ cascade by primary leptons with the monoenergetic (figures on the left) and the power law spectrum 
(on the right). 



(8 > 90°) can be again converted into the next generation 
of e pairs in the absorption process in the MSWR (see 
the difference between shock and esc in Table 1). Therefore, 
possible contribution of the cascades initiated by these sec- 
ondary e pairs in the MSWR to the total 7-ray spectrum 
escaping from the binary system can not be neglected. These 
secondary e pairs propagate in the MSWR along the paths 
which are determined by the structure of the magnetic field 
in the wind of the massive star. That is why, the next gener- 
ation of 7-rays produced by them move in directions which 
can completely differ from the initial angular distribution 
of 7-rays escaping from the PWZ. This more complicated 
cascade with the re-distribution of charged particles with 
respect to their initial directions is discussed in details in 
the next section. 

3.2 Gamma-ray light curves and phase resolved 
spectra 

In order to investigate in more details the angular depen- 
dence of the 7-ray emission produced in the cascade inside 
the PWZ, we calculate the 7-ray light curves which should 
be observed at different inclination angles of the binary sys- 
tem. The results are shown for the two, HE and VHE, energy 
ranges (i.e. 0.1—10 GeV and 10— 10 3 GeV) in the case of pri- 
mary leptons injected into the PWZ with the monoenergetic 
and power law distributions (see Figs. 4). The 7-ray lumi- 
nosities are calculated after normalization of the power in 
primary leptons to the rotational energy lost by the pulsar, 
L e — Lrot- Clear anticorrelation in the 7-ray curves is ob- 
served between these two energy ranges for both initial spec- 
tra of primary leptons. The VHE emission is mainly limited 
to the phases when the pulsar is in front of the massive star 



(<p = 0°). In contrast, the HE emission is much more uniform 
with significant decrease at phases around 0° at which the 
VHE emission is the strongest. The disappearance of 7-ray 
emission for large inclination angles of i > 90° — a (where 
a = 26.4° is the angular extend of the massive companion 
observed from the distance of the binary separation) and at 
the phase close to 180°, is connected with the total eclipse of 
the pulsar by the massive companion. The 7-ray spectra for 
selected phases of the pulsar on its orbit and the inclination 
angles mentioned in Figs. 4 are shown in Figs. 5. The 7-ray 
spectral shapes do not differ significantly below a few GeV 
(spectral index close to -1.5) for half of the pulsar phases, 
independently on the inclination angle of the binary system. 
The level of this emission vary only by a factor of 2-3. In con- 
trast, the VHE emission is limited to a relatively small range 
of phases with clear dependence on the inclination angle of 
the binary. 



4 MASSIVE STAR WIND REGION 

As we have shown above, a significant part of 7-rays from 
the PWZ which passed the termination shock is effectively 
absorbed in the MSWR by dense radiation field of the mas- 
sive star (see Table. 1). However, the next generation of e ± 
pairs is forced to follow the local magnetic field lines and 
their directions of propagation can change significantly with 
respect to the directions of their parent 7-rays. Therefore, 
the development of the cascade in the MSWR becomes much 
more complicated, determined by the complex magnetic field 
structure and the cascade is not further one-dimensional. 
The next generation 7-rays are usually produced at com- 
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pletely different angles than the initial directions of the PWZ 
cascade 7-rays. 

It is assumed in these cascade calculations that syn- 
chrotron energy losses of the cascading e pairs can be ne- 
glected with respect to their energy losses on the ICS. In 
fact, the comparison of energy densities of the magnetic field 
~ 2.5 x 10 16 eV cm -3 , for the surface magnetic field of the 
star B s = 10 3 G, and the thermal radiation ~ 8.5 x 10 17 eV 
cm -3 , for the temperature of the star T B — 1.36 x 10 5 K, 
show that the IC losses should clearly dominate over the syn- 
chrotron process. Further from the star the energy density 
of radiation drops with the square of the distance and the 
energy density of magnetic field drops with the forth power 
of the distance (see Eq. (17), note that the Alfven radius is 
very close to the stellar surface for the considered here pa- 
rameters) . Therefore, if the IC losses dominate at the stellar 
surface it have to dominate everywhere above the star. It 
has been shown in Bednarek (1997, see Eq. 4 in that pa- 
per) that the Bremsstrahlung energy losses of leptons in the 
wind of the massive star are negligible with respect to their 
IC losses in the Thompson (T) regime for the parameters 
of the wind and the star considered here and the Lorentz 
factors of leptons above the limiting the value 10 3 consid- 
ered in the cascade process. However in the present paper 
we also consider leptons with energies up to 10 6 MeV. It is 
necessary to discuss in a more detail the case of the IC losses 
in the Klein-Nishina (KN) regime. The IC losses of leptons 
with the Lorentz factors 7 in the T regime and the thermal 
radiation with temperature T s can be approximated by 



(f) 



1.3 x 1(T 22 7 2 T S 4 



MeV s" 



(12) 



We can estimate energy losses in the KN regime by putting 
into the above formula the Lorentz factor of leptons between 
the KN and T regimes, i.e. 7 « 2 x 10 4 for T s = 10 5 K. Then 
we get 



(-) 

V dt J KN 



5 x l(f MeV s" 



(13) 



Note that the logarithmic dependence of energy losses in 
the KN regime has been neglected. On the other hand, the 
bremsstrahlung energy losses can be approximated by 



(—) 

\ dt ) b r. 



1.4 x 10~ Lb NE MeV s" 



(14) 



where N is the density of particles in cm~ 3 and E is energy of 
leptons in MeV. For the extreme energies of injected leptons 
considered in the cascade process (E = W 6 MeV), we obtain 



(-) 

V dt J brcm 



1.4 x 10" iu 7V MeV s" 



(15) 



iFxom the comparison of (dE / dt)nN with (dE/dt)^, 
obtain the critical density of the matter N a 

Let's estimate a typical density of the massive star wind 
in the MSWR close to the stellar surface, 

N w = M/4ttR 2 V w , (16) 
which for the wind velocity V w = 10 8 cm s _1 and the 



mass loss rate M = 4 x 10 5 M yr \ gives ~8x 10 13 
cm -3 . So then, density of the massive star wind is 2-3 or- 
ders of magnitude lower than the critical density, iV cr , above 
which the Bremsstrahlung energy losses dominate over the 
IC losses in the KN regime. Therefore, we conclude that the 
Bremsstrahlung process can be safely neglected with respect 
to the IC process for the model considered in this paper. 

Below we describe the procedure applied for the cal- 
culation of the cascade processes with the presence of the 
magnetic field in the MSWR. 



4.1 Magnetic field structure inside MSWR 

The magnetic field in the wind of the massive star can have 
complicated structure. In the region very close to the mas- 
sive star surface it is characterized by dipolar component. 
At a certain distance, the radial component starts to domi- 
nate due to the presence of the ionized plasma, and at larger 
distances the magnetic field becomes toroidal due to the ro- 
tation of the massive star. The strength of the magnetic 
field as a function of distance from the center of the massive 
star can be described by the following equations (Usov & 
Melrose, 1992), 



B(r) ^B s x 



R s <r< R A , 
R A <r < R toI , (17) 



(«rot/«oo) (R s /(RAr)), Rtor < r, 



where R tor is the radius of the toroidal field defined by the 
rotation velocity of the massive star w rot and the star wind 
velocity at the infinity Voo, Rtor = -Rs«oc /^rot • In the case of 
the typical parameters of considered here for the WR star, 
the surface magnetic field is of the order of B B = 10 3 G 
and the v TO t ~ 0.1i?oo. The Alfven radius is located then 
at Ra = 1.12 x R s , and R toI = 10 x i? s . Therefore, the 
main volume of the MSWR, in which IC e pair cascade 
can efficiently develop, is dominated by the magnetic field 
with the radial structure. 



4.2 Propagation of leptons in MSWR 

The equations of motion of a charged particle in a constant 
uniform magnetic field (see Jackson 1962) are described in 
the coordinate system in which the vector of the magnetic 
field, B, is parallel to the axis OZ. Then, the coordinate of 
the charged particle and its velocity is given by, 

r(t) — ro + v«te3 + acosu>Bte2 + asinuiBtei, (18) 
v(t) — vntes + ujB(i{cosu>Btei + sinwBtei), (19) 

where e*3 is a unit vector parallel to the magnetic field line, 
e*i and $2 are the other orthogonal unit vectors, lob = 
eB/^mc — ecB/E is the Larmor frequency and ujbcl = v± 
is the perpendicular velocity of the charged particle with 
respect to the magnetic field line while vu denotes compo- 
nent along the field, and a is the Larmor radius. Using these 
equations we follow the charged particle by determining its 
position and velocity vector with respect to the local mag- 
netic field line. 

The step distance method is applied to follow the path 
of the particle in the magnetic field. The dimension of the 
step is chosen to be less than 0.01 of the Larmor radius in 
the local magnetic field. The optical depth for lepton on 
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the IC scattering is calculated by summing up contributions 
from every step up to the moment of fulfilling the condi- 
tion for simulating the place of interaction (given by Eq. 
9). Note that this procedure is much more complicated and 
time consumming than applied in the PWZ since now we 
have to determine the direction of lepton at every step. We 
calculate the energy of produced secondary 7-ray photons 
applying the procedure described in Sect. 3.1 (see Eq. 11). 
The energy of surviving leptons is equal to the difference 
between its initial energy and the energy of produced 7-ray 
photons. Therefore, after every scattering the conditions for 
propagation of leptons in the magnetic field change due to 
the change of their Larmor radius. The particle paths are 
followed up to the moment when leptons: (1) are cooled be- 
low 500 MeV, (2) move outside the sphere with the radius 
of 10 massive star radii (since the efficiency of IC scattering 
at larger distances becomes low), (3) fall onto the massive 
star surface, or (4) move back through the shock plane. The 
energies and directions of produced 7-rays are stored and 
used for analysis of their spectral and angular properties. 

4.3 Angular distribution of escaping gamma-rays 

Let us investigate the angular distribution of secondary 7- 
rays produced by secondary pairs which in turn orig- 
inate in the absorption process of 7-rays from PWZ cas- 
cade. These e pairs propagate in the magnetic field and 
anisotropic soft radiation. As before, we consider two cases 
of the monoenergetic and power law injection of primary lep- 
tons. In Figs. 6 and 7 the distribution of secondary 7-rays 
on the sky are shown in the case of primary leptons injected 
at specific directions defined by the angles 9 = 90° (Figs. 6a 
and 7a), 120°(Figs. 6band7b), and 150° (Figs. 6c and 7c). 
These numbers are shown for one injected monoenergetic 
lepton with energy 10 6 MeV (Fig. 6) and for the power law 
spectrum of leptons with the power normalized to 10 6 MeV 
(Fig. 7). The 7-rays are sorted by their escape directions, i.e. 
within the regions on the sky defined by the azimuthal di- 
rections (ip = 180° corresponding to the case of the massive 
star in front of the pulsar) and by the cosine of the vertical 
angle j3. The basic features of the angular distribution of 
the secondary 7-rays are determined by the structure of the 
magnetic field in the MSWR. For the parameters considered 
in this work, the magnetic field is mainly radial in the re- 
gion above the shock. For the primary leptons injected at the 
angle 9 = 90°, the 7-rays from cascade inside the MSWR 
are mainly produced inside the cone centred along the di- 
rection of the magnetic field lines which cross the direction 
of propagation of 7-rays produced inside the PWZ. Since 
the magnetic field lines are radial, the angles at which the 
secondary 7-rays escape are limited to relatively small part 
of the sphere (see Figs. 6a and 7a). For primary leptons in- 
jected at larger angles, see e.g. 9 — 120° and 150° (Figs. 6bc 
and 7bc), the cone of secondary 7-ray production inside the 
MSWR becomes broader on the sky. Most efficient produc- 
tion occurs for directions which are tangent to the massive 
star limb (see Figs. 6c and 7c for the angle 9 — 150°). In 
this case a significant part of secondary 7-rays fall onto the 
surface of the massive star which angular dimensions, seen 
from the location of the neutron star, are equal to 26.4° (note 
the avoidance region centred on tp = 180° and cos/3 = in 
Figs. 6c and 7c). If the primary leptons are injected exactly 




1 h I 




<P 

Figure 6. The maps with the numbers of 7-ray photons with en- 
ergies > 100 MeV (in log grey scale) escaping at specific directions 
on the sky defined by the azimuthal angle <p and the cosine of the 
vertical angle (i (the bins with A(p = 10° and cos fi = 0.1). ip = 0° 
correspond to the case of the massive star in front of the pulsar 
which is the source of primary leptons. These 7-rays are produced 
in the cascade in the MSWR initiated by the secondary cascade 
7-rays escaping from the PWZ, which in turn were produced by 
primary monoenergetic leptons with energies 10 6 MeV. The re- 
sults, averaged over N simulated primary leptons, are shown for 
different initial injection angles of the primary leptons equal to 
6 = 90° for N = 200 (a), 120° for N = 100 (b), and 150° for 
N = 60 (c), which are measured above the plane of the binary 
system (i.e. towards positive vertical angles). Presented results 
are re-normalized for one lepton with energy 10 6 MeV. 



towards the massive star (9 = 180°), then the secondary e 
leptons, which originate in the MSWR in absorption of sec- 
ondary 7-rays from PWZ, move along the radial magnetic 
field. Most of the 7-rays produced in MSWR fall then onto 
the massive star. The secondary e pairs reach finally the 
dipole part of the magnetic field, which is close to the mas- 
sive star, and may produce low energy 7-rays at wide angles 
in the plane of the binary system. 

The angular distribution of 7-rays produced in the 
MSWR in the case of monoenergetic and power law distribu- 
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Figure 7. As in Fig. (6) but for the primary leptons with the 
power law spectrum (described in the text) which energy is also 
normalized to 10 6 MeV. The results are averaged over 500 simu- 
lated primary leptons for the angle of injection 9 = 90° (a), 180 
leptons for 120° (b), and 150 leptons for 150° (c). 



tions of primary leptons is similar since it is determined by 
the radial structure of the magnetic field. However, the num- 
bers of produced 7-rays are a factor of a few larger in the case 
of monoenergetic injection of primary leptons. This is due 
to the fact that 7-rays produced in the PWZ by monoener- 
getic leptons have on average higher energies (see Figs. 3ac). 
Therefore they transport more energy to the MSWR. The 
energies and numbers of secondary pairs, from their ab- 
sorption inside the MSWR, are larger allowing more efficient 
production of next generation of 7-rays in the MSWR. 

The angular distributions of 7-rays with energies > 100 
MeV which escape to the observer from the MSWR for the 
case of isotropic injection of primary leptons are shown in 
Fig. 8. The numbers of these 7-rays are normalized to the 
power in primary lepton spectrum equal to 10 6 MeV sr _1 . 
As before, the pulsar is behind the massive star for the phase 
if — 180°. As expected, the numbers of produced 7-rays are 
the largest in directions tangent to the limb of the massive 
star. Significantly lower numbers of 7-rays emerge from the 
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Figure 8. The map with the total numbers of 7-ray photons 
with energies > 100 MeV (in linear grey scale), escaping from 
the binary at specific directions on the sky (defined by the az- 
imuthal angle <p and cosine of horizontal angle cos 8) for the case 
of isotropic injection of primary monoenergetic leptons with en- 
ergies 10 6 MeV (a) and for the primary leptons with the power 
law spectrum described in the text (b). 7-rays are produced in 
cascades inside the MSWR initiated by secondary cascade 7-rays 
from the PWZ. The spectra of primary leptons are normalized to 
their total energy equal to 10 6 MeV per steradian. 



binary system in the outward directions, i.e. in directions for 
which the pulsar is in front of the massive star with respect 
to the observer. Note that in these directions the MSWR 
does not extend (see Fig. 1). 



4.4 Gamma-ray light curves and spectra 

The 7-ray light curves of photons produced in the MSWR 
(Table 4) for two energy ranges, 0.1 — 10 GeV and 10 GeV 
- 1 TeV, are shown in Figs. 9, after its normalization to 
the total energy loss rate of the considered here pulsar. The 
observer is located at different inclination angles i to the 
plane of the binary system. The general features of the 7- 
ray light curves produced in the MSWR are similar to these 
ones shown above for the 7-rays produced in the PWZ (see 
for comparison Figs. 9). However, the total 7-ray power is 
usually lower (typically by a factor of three) than expected 
in the case of 7-rays produced in the PWZ. Only in the 
range of phases around 0°, and large inclination angles, the 
power of the 7-rays from the MSWR can dominate over the 
PWZ 7-rays in the GeV energy range. 

The power emitted in GeV energies is quite uniform 
with the phase of the pulsar provided that the inclination 
angle is small enough that the binary system is not eclipsing. 
The power emitted at TeV energies is limited to the range of 
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log(E / MeV) log(E / MeV) 

Figure 10. The 7-ray spectra escaping from the MSWR for different phases of the pulsar with respect to the location of the observer, 
ip = 30° (solid curve) 60° (dashed), 90° (dot-dashed), 120° (dotted), 180° (dot-dot-dot-dashed), and three ranges of the cosine of the 
inclination angle of the binary 0.76 < cos a < 0.88 (centered on i = 30° - (a) and (d)), 0.41 < cos a < 0.53 (60° - (b) and (e)), and 
—0.06 < cos a < 0.06 (90° - (c) and (f)). 7-rays are produced in the MSWR cascade by secondary e ± pairs which in turn originate in 
the absorption of 7-rays produced in the PWZ by primary leptons with the monoenergetic (figures on the left) and power law spectrum 
(on the right). 




60 120 180 240 300 60 120 180 240 300 360 

[•] <t> ['] 

Figure 9. The 7-ray light curves in two energy ranges, 0.1 — 10 
GeV ((a) and (c)) and 10 - 10 3 GcV ((b) and (d)), produced in 
the MSWR cascade initiated by secondary 7-rays which are in 
turn produced in the cascade in the PWZ initiated by primary 
leptons with the monoenergetic (figures on the left) and the power 
law distributions (on the right). 7-rays escaping from the MSWR 
are collected within the range of the cosine of inclination angles 
of the system 0.76 < cos i < 0.88 (contains i = 30° - solid curve), 
0.41 < cosi < 0.53 (60° - dashed), and -0.06 < cosi < 0.06 (90° 
- dotted). 



phases ±60° around the phase 0° which correspond to the 
position of the pulsar in front of the massive star. 

Figs. 10 show the 7-ray spectra produced in the MSWR 
for selected phases and inclination angles of the binary sys- 
tem. The basic differences between these spectra and the 



7-ray spectra produced in the PWZ are due to the fact that 
cascade in the MSWR is initiated by leptons for which the 
total optical depths are higher. Therefore, these 7-ray spec- 
tra are steeper, with the spectral index close to -2, in con- 
trast to the spectral index of the spectra produced in the 
PWZ which are close to -1.5. Although, the power emitted 
in 7-rays from the MSWR is a factor of a few lower. It can 
still significantly contribute to the total 7-ray spectrum ob- 
served at some range of phases and at energies below ~ 1 
GeV for the case of monoenergetic injection of primary lep- 
tons. Moreover, as already noted above, this GeV emission 
is much more uniform over the sky. At energies > 10 GeV 
the contribution of the 7-rays from the MSWR (produced by 
monoenergetic primary leptons) becomes negligible with re- 
spect to the 7-rays from the PWZ for the considered location 
of the shock inside the binary system defined by r\ — 0.3. The 
7-rays produced inside the MSWR by primary leptons with 
the power law spectrum (see Figs. lOdef) do not contribute 
significantly to the 7-rays produced in the PWZ (Figs. 5def). 
Primary leptons with the power law spectrum produce on 
average more secondary 7-rays already inside the PWZ. 



5 DEPENDENCE ON THE SHOCK 
LOCALIZATION 

Since detailed parameters of the massive stellar winds are 
not precisely known we investigate the dependence of es- 
caping 7-ray spectra on the location of the shock inside the 
binary system by performing calculations for the same initial 
injection spectra of leptons but for the shock determined by 
the value of ij = 0.06. This value is obtained for upper limits 



© RAS, MNRAS 000, 000-000 



Gamma-rays from cascades in close massive binaries containing energetic pulsars 13 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

- (a) ■ ' (c) 




60 120 180 240 300 360 60 120 180 240 300 360 



<t> ['] * f] 

Figure 11. As in Fig. 4 but for r? = 0.06. 

of the stellar wind parameters (see Sect. 2). For 77 = 0.06 the 
shock is located at the distance of p = 0.25 i? s from the 
pulsar, i.e. about two times closer than for ij = 0.3, and the 
volume for the development of the cascade in the PWZ is 
much smaller than for the cascade in the MSWR. 



5.1 Gamma-rays from PWZ 

For leptons injected with the power law spectrum, the 7-ray 
luminosities produced in the PWZ (and detected by the ob- 
server at the infinity) for 77 = 0.06 are only slightly lower 
than for 77 = 0.3 (see for comparison the 7-ray light curves 
and spectra in Figs. 11 and 12 with Figs. 4 and 5, respec- 
tively). This small differences are not surprising since the 
energy from primary leptons with the power law spectrum 
is transfered to 7-rays efficiently already inside the PWZ 
due to large optical depths for leptons with energies much 
lower than 10 6 MeV (see Fig. 2). However, the reduction of 
7-ray fluxes is much stronger for leptons injected with the 
monoenergetic spectrum and energies 10 6 MeV. The 7-ray 
light curves and spectra for the monoenergetic primary lep- 
tons and the shock defined by 77 = 0.06 are by a factor of 
about 3 lower than for 77 = 0.3. As shown in Fig. 2c, the 
optical depths for leptons with energies 10 6 MeV are close 
to unity for 77 — 0.06 but much larger for 77 = 0.3. Therefore, 
the differences in the cascades developed in the PWZ by 
high energy monoenergetic leptons are larger than reported 
for leptons with the power law spectrum. 

Large differences between these two discussed localiza- 
tions of the shock in the 7-ray light curves in the energy 
range 10 — 10 GeV are also observed for small inclination 
angles i (see Figs. 4b and lib). They are due to lower opti- 
cal depths for monoenergetic leptons with energies 10 6 MeV. 
Leptons injected in directions perpendicular to the plane of 
the binary system meet on its propagation through the PWZ 
smaller number of soft photons from the massive companion 
than these ones propagating towards the massive star due to 
the geometrical effects. Moreover, the spectra of 7-rays es- 
caping to the infinity for 77 = 0.06 show cut-offs at lower en- 
ergies than in the case of the shock defined by 77 = 0.3. This 
effect, especially clear for the monoenergetic injection of pri- 
mary leptons, is due to the stronger absorption of higher 



energy 7-rays during their propagation in the MSWR which 
volume is larger for 77 = 0.06 than for 77 = 0.3. 

5.2 Gamma-rays from MSWR 

Basic differences in the 7-ray production in the MSWR due 
to different localizations of the shock inside the massive bi- 
nary (defined by 77 = 0.3 and 0.06) are related to the ef- 
ficiency of the cascades occurring already inside the PWZ 
since e ± pairs, created in absorption of the PWZ 7-rays, 
are responsible for production of 7-rays in the MSWR. For 
localization of the shock described by 77 = 0.06, the 7-ray 
spectra produced in the MSWR extend to higher energies 
than for the case of 77 = 0.3 (see Figs. 10 and 14). This is 
also clearly seen when investigating the 7-ray light curves 
at energies > 10 GeV for the case of the power law spec- 
trum of primary leptons shown in Figs. 9d and 13d. This 
behavior is in contrast to the features of the 7-ray spectra 
produced in the PWZ. It is clear from the Table lb that, al- 
though the 7-ray luminosity leaving the PWZ (i.e. entering 
the MSWR) for the power law spectra of primary leptons 
is comparable for 77 = 0.3 and 0.06 (see percentages marked 
by shock), the power transformed to the secondary e pairs 
in the MSWR, from absorption of the PWZ 7-rays, differ 
significantly due to much larger volume of the MSWR for 
77 = 0.06 (see differences between percentages marked by 
shock and esc). Therefore, secondary e pairs produced in 
the MSWR for the shock located closer to the pulsar (i.e. 
for 77 = 0.06) have, on average, higher energies and produce 
more energetic 7-rays in the MSWR. This is also the reason 
that the 7-ray light curves are on a slightly higher level for 
77 = 0.06 (compare Figs. 9cd with 13cd). 

For the case of the monoenergetic injection of primary 
leptons, the 7-ray luminosity originating in the PWZ is 
clearly lower for 77 = 0.06 (Table la). However, the energy 
transfered into the secondary e pairs in the MSWR, as a 
result of absorption of the PWZ 7-rays, is comparable for 
both of the shock localizations (note differences between per- 
centages marked by shock and esc in Table la). Therefore, 
differences in the 7-ray light curves and spectra for the lo- 
calizations of the shock inside the binary system defined by 
77 = 0.3 and 0.06 are not so evident than in the case of the 
power law injection of primary leptons. 



6 DISCUSSION AND CONCLUSIONS 

We have analyzed in details the anisotropic IC e ± cascade in- 
side the massive binary system containing young pulsar. We 
have considered the situation in which leptons are injected 
mono-directionally or isotropically from a place (identified 
with the pulsar) which differs from the center of the isotropic 
source of soft radiation (identified with the massive stellar 
companion). The simplest geometrical situation has been 
analyzed in which a compact object is on a circular orbit 
around the massive star. As an example, the parameters of 
the massive binary in the Cyg X-3 system have been applied. 
We concentrate on the situation in which a pulsar is ener- 
getic enough to create strong shock in collision with the wind 
of the massive star. As a consequence, two regions for the 
development of the cascade has to be distinguished which 
differ in basic properties, i.e. the pulsar wind zone (PWZ) 
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and the massive star wind region (MSWR). In the PWZ 
the cascade develops radially from the pulsar and therefore 
can be considered only in one direction. In the MSWR the 
cascade develops in complex magnetic field of the massive 
star. Therefore, complicated propagation paths of secondary 
e pairs in this magnetic field have to be taken into ac- 
count. However, for typical values of the magnetic fields in 
the winds of the massive stars and parameters of Cyg X- 
3 binary, the energy density of radiation from the massive 
star is larger than the energy density of the magnetic field. 
Hence, the IC energy losses of leptons dominate over their 
synchrotron losses and the influence of synchrotron process 



on the energy budget of leptons can be neglected in the first 
approximation. 

For this specific binary system we have calculated the 
optical depths for leptons and 7-rays and found that they 
are much larger than unity (Fig. 2). Therefore, if relativis- 
tic leptons are injected inside the binary they should ini- 
tiate efficient cascades. Such cascades are analyzed in the 
case of leptons injected by the compact object with: (a) the 
monoenergetic spectrum and energies 10 6 MeV and, (b) the 
power law spectrum of the type predicted by the calculations 
of the cascade processes in the inner pulsar magnetosphere 
(Hibschman & Arons 2001). 
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Figure 13. As in Fig. 6 but for r\ = 0.06. 

The main results of our calculations are following: 

(i) Most of the energy of primary leptons is transferee! to 
the secondary 7-rays (~ 80 — 90%, see Table 1) for both the 
monoenergetic and the power law primary spectra, and for 
the shock defined by r\ — 0.3 which correspond to the case of 
the massive star with the wind parameters from the middle 
of the estimated range (Sect. 2). This percentage of energy 
does not depend strongly on the location of the shock for the 
power law spectrum of primary leptons (compare ij — 0.3 
and 0.06 in Table 1), but is significantly lower for the shock 
closer to the pulsar, determined by r\ — 0.06 in the case 
of monoenergetic primary leptons. The difference between 
the power in 7-rays arriving to the shock (marked by shock) 
and 7-rays escaping to the observer (marked by esc) gives 
part of the power of 7-rays produced in the PWZ which is 
converted into the secondary pairs in the MSWR. These 
e pairs can initiate additional cascade in the MSWR and 
take typically between one third and half of the 7-ray power 
escaping to the observer from the PWZ (see Table 1). There- 
fore, the contribution of 7-rays produced in the MSWR is 
usually lower than 7-rays produced in the PWZ. However, 
if the monoenergetic primary leptons are injected into the 
binary with the shock located relatively close to the pul- 
sar (e.g. for rj — 0.06) the 7-ray luminosity produced in the 
MSWR can be even a factor of two or three larger than 7- 
ray luminosity from the PWZ. This is due to lower optical 
depths in relatively small PWZ and so inefficient cascading 
precess initiated by high energy monoenergetic leptons. 

(ii) The 7-ray light curves for photons escaping from the 
PWZ at lower energies (energy range 0.1 — 10 GeV) are an- 
ticorrelated with the 7-ray light curves at higher energies 
(> 10 GcV) for both discussed primary spectra of leptons 
independently on localizations of the shock within the bi- 
nary system (see e.g. the 7-ray light curves in Figs. 4 and 

11) . This basic feature is also clearly seen when comparing 
the 7 spectra escaping to the observer at different phases of 
the pulsar and inclinations of the binary system (Figs. 5 and 

12) . It is easily understood since the intense 7-ray fluxes at 
energy range 0.1 — 10 GeV are expected in directions of large 
optical depths, when the binary system is viewed from the 
direction close to the massive star. In contrast, high level 
7-ray fluxes at energies > 10 GeV are expected when the 



compact object is in front of the massive star. For consid- 
ered parameters of the binary system, the optical depth is 
high enough for efficient production of 7-rays above 10 GeV 
but is too low for efficient cascading and production of lower 
energy 7-rays. Moreover, the 7-ray spectra extend to higher 
energies for the monoenergetic primary leptons but are less 
intense at lower energies than expected for the power law pri- 
mary leptons (see Figs. 5 and 12) since the optical depths are 
lower for higher energy particles (see Figs. 2) . Therefore, it is 
expected that detailed results in the 7-ray spectra from spe- 
cific binaries (PSR B1259-63, SAXJ 0635+0533), for which 
the parameters of the massive star are well known, should 
allow to obtain information on the spectrum of primary lep- 
tons accelerated by the pulsar. This problem is difficult to 
investigate directly in the case of isolated pulsars for which 
the surrounding soft radiation field in the pulsar wind zone 
is relatively weak. 

(iii) The complex magnetic field in the MSWR has sig- 
nificant effect on the directions of propagation of secondary 
e pairs and, as a consequence, on the distribution of 7-rays 
on the sky. Our detailed calculations show that the angular 
distribution of 7-rays, produced in the cascade processes in 
the MSWR by secondary e pairs from absorption of 7-rays 
from PWZ, strongly depends on the injection direction of 
primary leptons (see maps in Figs. 6 and 7). For example, 
if the primary particles are injected perpendicular to the 
plane of the binary system then the 7-rays produced in the 
cascade are collimated within the cone around the direction 
of the local magnetic field lines, i.e. in the outward direc- 
tion with respect to the massive star (Figs. 6a and 7a) . This 
is due to the strong radial component of the massive star 
magnetic field in the main volume of the binary system (see 
Eq. 11). Therefore, in the case of highly anisotropic injection 
of primary particles by the pulsar, the strong 7-rays fluxes 
may appear in quite different directions. In contrast, if the 
primary particles are injected in the directions not far from 
the the stellar limb, this radial component of the magnetic 
field tends to focus the 7-rays in directions tangent to the 
stellar limb. A significant part of these 7-rays impinge onto 
the surface of the massive star (see Figs. 6c and 7c) . It is 
expected that these 7-rays can excite observable nuclear 7- 
ray line fluxes due to their interactions with the matter on 
the massive star surface. This line emission should be corre- 
lated with the high energy 7-ray emission (> 10 GeV) and 
anticorrelated with the low energy 7-ray emission (between 
0.1 - 10 GeV). 

The distribution of 7-rays produced in the MSWR is also 
strongly anisotropic for the case of isotropic injection of pri- 
mary leptons by the pulsar (see Fig. 8). There is a strong 
concentration of produced secondary 7-rays in directions 
close to the limb of the massive star. The number of sec- 
ondary photons produced by the primary leptons with the 
power law spectrum is by a factor of a few lower than for mo- 
noenergetic primary leptons but general emission pattern on 
the sky does not differ significantly. These features are due 
to the fact that leptons with lower energies produce 7-rays 
with energies above 100 MeV less efficiently. 

(iv) The 7-ray fluxes produced in the MSWR are weakly 
dependent on the phase of the pulsar and the inclination 
angles of the binary system, except for the case of eclips- 
ing binaries (see Figs. 9 and 13). This is due to the prop- 
agation effects of secondary e ± pairs in the magnetic field 
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of the MSWR which result in a significant isotropisation 
of the secondary e pairs with respect to the directions of 
propagation of primary leptons. Therefore, when the pulsar 
is in front of the massive star (phase ip = 0°), the main 
contribution to the escaping low energy 7-ray flux comes 
from the MSWR and the anticorrelation between the lower 
energy (0.1 — 10 GcV) and higher energy (> 10 GeV) 7- 
rays produced in the PWZ is partially suppressed (compare 
Figs. 4 with 9 and 11 with 13). The 7-ray spectra escap- 
ing to the observer from the MSWR are flatter in the GeV 
energy range (spectral index close to -2) than the 7-ray spec- 
tra from the PWZ (spectral index close to -1.5) (compare 
Figs. 5 with 10 and 12 with 14). This is due to the difference 
in the optical depth for particles propagating in the PWZ 
and the MSWR. Since the cascade is not efficient enough in 
the PWZ, the secondary e ± pairs are not able to influence 
lower energy part of the IC 7-ray spectrum formed mainly 
by cooling of primary leptons. However, secondary e pairs 
in the MSWR are cooled to energies low enough to produce 
a lot of secondary GeV photons and thus creating steeper 
spectra. Due to the same reasons, the 7-ray spectra from the 
MSWR extend usually to lower energies and the predicted 
fluxes above ~ 10 GeV are lower. 

(v) The 7-ray spectra escaping from the MSWR depend 
strongly also on the localization of the shock inside the bi- 
nary. This is caused by the change in relative volume for 
the development of cascades in the PWZ and MSWR. If the 
shock is closer to the pulsar then the 7-ray spectra from the 
MSWR extend to higher energies (see Figs. 10 and 14) since 
they are produced by on average more energetic secondary 
e ± pairs originated in the absorption of more energetic 7- 
rays escaping from the PWZ. The 7-ray fluxes at energies 
< 10 GeV, produced in the MSWR, do not depend strongly 
on the localization of the shock. But the 7-ray fluxes at en- 
ergies > 10 GeV are on a much lower level for the case of 
primary leptons injected with the power law spectrum (see 
Figs. 9 and 13) from the same reasons as discussed in this 
item above. 

In summary, detailed investigation of the phase resolved 
IC continuous 7-ray spectra and the 7-ray light curves from 
the compact massive binaries, in which energetic pulsar is re- 
sponsible for the acceleration of relativistic leptons, should 
allow extraction of information on the primary spectra of 
leptons (monoenergetic, power law ?) and conditions for col- 
lisions of the pulsar and stellar winds. The fluxes and spectra 
of escaping 7-rays should be also sensitive to the acceleration 
site of the primary leptons, inside the inner pulsar magne- 
tosphere (power law spectrum ?), close to the pulsar light 
cylinder (monoenergetic ?), or during the propagation in 
the pulsar wind zone (see e.g. re-acceleration model of Con- 
topoulos & Kazanas 2002) . The effects on the 7-ray spectra 
caused by leptons which are additionally accelerated in the 
pulsar wind zone have been recently investigated with the 
application to the PSR B1259-63/Be binary system (Sier- 
powska & Bednarek 2004b), but they need further more de- 
tailed calculations. 

Taking into account sensitivities of the future and al- 
ready operating 7-ray telescopes in the energy range below 
a few hundred GeV (which are estimated on ~ 3 x 10~ 13 erg 
cm" 2 s" 1 for the GLAST and ~ 10~ 12 erg cm" 2 s" 1 for the 
MAGIC), the type of the massive binary considered in this 



paper with the 7-ray luminosity of the order of ~ 10 erg 
sr _1 s _1 (similar to the Cyg X-3 massive binary) should be 
detected at any location within our Galaxy and also within 
the Magellanic Clouds. 

In the forthcoming paper (Sierpowska & Bednarek 
2004a) we intend to apply such general model for the 7-ray 
production in specific massive binaries in which case there 
are evidences of the existence of energetic pulsars (e.g. PSR 
B1259-63, SAX J0635) or compact objects with jets (see 
the list of sources in Table 1 in Mirabel & Rodriguez 1999). 
Note that in the case of less compact binaries than con- 
sidered here, the secondary cascade e ± pairs arrive to the 
shock with significant energies since the optical depths for 
leptons are much lower. These pairs move then along the 
shock plane. Their contribution to the total escaping 7-ray 
fluxes has to be also taken into account. It has not been con- 
sidered in this paper since, for the parameters of the massive 
star in the Cyg X-3, these pairs take only a relatively small 
part of energy of primary leptons (typically less than ~ 20% 
of the total injected energy, see Table 1). Therefore, relative 
importance of processes occurring in the PWZ, the shock 
region, and the MSWR can differ significantly for the very 
close binaries with the pulsars on circular orbits (as consid- 
ered here) from the binaries of the PSR 1259-63 type, i.e. 
with a broad, highly eccentric orbit and a massive Be star 
producing non-spherical stellar wind. 

We also plan to investigate possible importance of the 
cascade processes inside the binary systems of two very mas- 
sive stars which create strong shock as a result of their stellar 
wind collisions, e.g. 7-2 Vel (WC8+07.5). Energetic e^ 1 pairs 
can appear inside such binaries as a result of interaction of 
hadrons, accelerated by the shock, with radiation field of the 
massive star. 
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